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S Abstract: We consider a quantum particle coupled (with strength A) to a spatial array of independent non-interacting reservoirs in thermal 
states (heat baths). Under the assumption that the reservoir correlations decay exponentially in time, we prove that the motion of the particle 
is diffusive at large times for small, but finite A. Our proof relies on an expansion around the kinetic scaling limit (A \ 0, while time and space 
scale as A~ 2 ) in which the particle satisfies a Boltzmann equation. We also show an equipartition theorem: the distribution of the kinetic energy 
of the particle tends to a Maxwell-Boltzmann distribution, up to a correction of 0(A 2 ). 
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2 ■ 1 Introduction 

oc : 

O ■ l.i Diffusion 



Diffusion and Brownian motion are among the most fundamental phenomena described by transport theory. They 
refer to the apparent random motion of a particle or, for that matter, any degree of freedom, interacting with 
many other, mutually independent degrees of freedom in a thermal state. The interactions produce an erratic 
macroscopic motion that we perceive as diffusive or as Brownian motion. From a mathematical point of view, we 
may attempt to understand diffusive motion by invoking a central limit theorem: N interactions produce an effect 
Sx, which is given by Sx ~ y/N. Since the number of interactions is proportional to the time lapse St, we can write 
(Sx) 2 = DSt, where the proportionality constant D is called the diffusion constant. Via the Einstein relation, the 
diffusion constant determines quantities such as the thermal or electric conductivity. 

The model of a particle (quantum or classical) coupled to a thermal reservoir of free particles is a natural starting 
point for an analysis of diffusion. We assume the particle to be quantum mechanical. By ( )^ we denote the 
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expectation value in a state where the reservoir has an inverse temperature < oo. Then 

({8xf)p = ((x(t) - x(0)) 2 ) p = ( dh [ dt 2 (x(h)x(t 2 ))0, (1.1) 

Jo Jo 

where xlt) is the position of the particle at time t and x(t) — i[H, x{t)], where H is the Hamiltonian of the system, 
is the velocity. We expect that, because of interactions with the reservoir, x(ti) and £(£2) become de-correlated 
rapidly, as |<i — £2! grows. Thus, the quantity | (±(ti)i(t2)) p\ is expected to be integrable in the variable ti —t\. 
Combining this with isotropy, i.e., (x(t))p —> rapidly, as t — ► 00, for (3 < 00, one concludes that, asymptotically as 
t tends to 00, 

((x(t)-x(0)) 2 ) p = D({3)\tl (1.2) 
for some positive, finite constant < D((3) < 00, given by 

D((3) = f dt(x(0)x(t)) p . (1.3) 

JR 

Because of the equipartition theorem, one expects that 

{x(t) 2 ) p ^ f dvv 2 e~^ v \ ast/oo, (1.4) 

where E(v) is the kinetic energy of a particle with velocity v. Obviously, i 1.4| | is strictly positive for finite (3. 
Likewise, we expect that D(/3) is strictly positive, for (3 < 00. 

Equations (|1.2Jl and (|1.4) l suggest that, at very large times, the motion of a particle interacting with a reservoir 
or heat bath at strictly positive temperature has universal features: The mean value of its speed is strictly positive 
and finite, and its mean displacement is proportional to the square root of time. In contrast, at zero temperature 
(13 = 00), the nature of the particle's motion depends on properties of the reservoir and the dispersion law, e{k), of 
the particle; (k £ is its momentum). If, for a particle momentum k, 

e(k- q) +uj{q) > e(k), forallg^O, (1.5) 

where is the dispersion law of a mode (particle) of the reservoir with momentum q, then the particle cannot 
lower its energy and reduce its speed by exciting a reservoir mode, i.e., by spontaneously emitting a reservoir 
particle. Its motion will therefore be ballistic. The only effect of the reservoir is a renormalization of the effective 
mass (the dispersion law e) of the particle. If, however, i 1.51 1 is not satisfied, then the particle can excite reservoir 
modes (emit reservoir particles). This process reduces its kinetic energy and speed, i.e., it leads to friction. Friction 
takes place at all momenta k if, e.g., Ld(q) oc \q\ 2 (reservoir particles are non-relativistic). If io(q) = c\q\, i.e., the 
reservoir particles are low-energetic phonons or photons, friction only takes place at momenta k of the particle 
where | Ve(fc)| > c. The radiation corresponding to the reservoir particles emitted in the process of friction is called 
Cerenkov radiation. 

Despite the importance of diffusion and its conceptual simplicity, there has, so far, not existed any rigorous 
proof that it occurs in a model as described above. In the present paper, we establish diffusion for models where 
the particle is coupled to a spatial array of independent heat baths. 

1.2 Informal description of the model and main results 

We consider a quantum particle hopping on the lattice Z d . With each lattice point, we associate an independent 
thermal reservoir consisting of a free bosonic quantum field describing phonons or photons at temperature 
(In this section, we present a description of the system appropriate at zero temperature; it is formal when (3 < 00.) 
The total Hilbert space, of the coupled system is a tensor product of the system space, J%$, with a reservoir 
space, J^r, which is a (separable) subspace of the infinite tensor product of reservoir spaces , x £ Z d , at all 
sites. Thus 

:= J^ s ® Jdk- (1-6) 



2 



The system space J£§ is given by l 2 (Z d ), and the particle Hamiltonian is given by the finite-difference Laplacian 
A. Each reservoir is described by a boson field; creation and annihilation operators creating /annihilating bosons 
with momentum q e M. d at site x are written as a x (q),a x (q) respectively, and satisfy the canonical commutation 
relations 

[a*(q),a*tf)] = 0, [a x (q), a* x , (</)] = 6 x , x >S(q - q'), (1.7) 

where a# stands for either a or a* . 

The total Hamiltonian of the system is taken to be 

H X :=-A+J2 I dqw(q)al(q)a x (q) + xY] [ dq \x)(x\ <8> {#?)(£ (?) + h.c.} , (1.8) 

where <f)(q) is a form factor and A e R is the coupling strength. We are writing A instead of A® 1 and a x (q) instead 
of 1 ® a x (q) 

The independence of the reservoirs has far-reaching consequences. Consider the lattice translation T z , z S Z d , 
acting on operators on by 

T z (\x)(y\) := \x + z){y + z\ (1.9) 
T,(o#(g)) := at + M- d-10) 

It is easily seen that 

T Z (H X ) = H X . (1.11) 

Notice that this transformation does not involve the momentum coordinates q inside the reservoirs. It is the exis- 
tence of this translation symmetry that allows us to obtain results on diffusion without very hard work. Assume 
we had started from a model with only one reservoir, with Hamiltonian given by 

H x :=-A+ [ dqcj{q)a*(q)a(q) + A V / dq \x){x\ ® {a^^^e" 1 ^ + h.c.) , (1.12) 

X £Z d 

where, now, the operators a(q), a*(q) do not carry an index x and (•, •) is the scalar product on C d . This model 
still exhibits translation symmetry, but this symmetry maps a*(q) — * a* (q)e~^ z ' q \ a(q) — > e^ z ' qS> a(q), which is 
the reason for the factor e~ l ( x -i) in the interaction Hamiltonian of ( |1.12t and leads to bad decay properties of the 
reservoir correlation functions. 

The initial state for the reservoirs is chosen to be '■= ® x ez d PR > where each is an equilibrium state at 
inverse temperature (3 for the reservoir at site x. For mathematical details on the construction of infinite reservoirs, 
see 1 7, 3, 1 ]. In Lemma l2~3l we define the reduced Heisenberg-picture dynamics (i.e., the particle dynamics obtained 
by tracing out the reservoir degrees of freedom) 

S i ► Z^*{S) := p f i [e UH *{S ® l)e- 1 *^] , S e &(Jt%) (1.13) 

Placing the particle initially at site 0, that is, in the vector |0), we study the distribution function 

rf(x) := (Q\Z^(\x){x\)\Q). (1.14) 

The quantity fJ-t(x) > is the probability to find the particle at site x at time t. One easily checks that 

E Vtix) = 1, (1-15) 

xGZ d 

and hence it is justified to think of fx* (•) as a probability density on Z d . By diffusion, we mean that, for large t, 

rf(x) - (^) d/ "(detC A )- 1 /2 cxp{ _ ^_L )jD -i^ }) ( L16 ) 
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where D\ = D\{j3) is a positive-definite matrix with the interpretation of a diffusion tensor. (Actually, if the 
particle Hamiltonian is given by —A (as in this section), the tensor D\ is isotropic and hence a scalar). We now 
move towards quantifying (1.16) . Let us fix a time t. Since Ht( x ) is a probability measure, one can think of xt as a 
random variable such that 

Prob A (xt = x) :=v$(x). (1.17) 

The claim that the random variable ^ converges in distribution, as t / oo, to a Gaussian random variable with 

mean and variance D^ 1 is called a Central Limit Theorem (CLT). It is equivalent to pointwise convergence of the 
characteristic function, i.e., 



V e~^ {x ' q) ^(x) — ► e -v(a,Dx9) ) forallg e R d , (1.18) 

c6Z d 



t]oo 

and it is this statement which is our main result, Theorem |3.2| 

Let X := J2xei, d x \ x )( x \ ^ e tne position operator on the lattice and write X t := Z^'*(X). Then a slightly 
stronger version of dl-181 1 implies that 

(0|^|0) — > 0, (0|^|0) — > D A , (1.19) 

and this will also follow from our results; see Remark l3.3[ 

Our second result concerns the asymptotic expectation value of the kinetic energy of the particle. Let E t := 
Z^'* (—A) be the kinetic energy at time t. We prove that, for all bounded functions 9, 

<» iw)i »> s Ji ix-"-~r' +°<^ 

where e(fc) = 5^,- =1 (2 — 2 cos fc J ) is the dispersion law of the particle. This is stated in Theorem |3.1| 



1.3 Related results 

In the physics literature, the model with Hamiltonian (1. 12) and with reservoir particles being phonons is referred 
to as the polaron model. We refer to E31I2D and references therein for a discussion. The first rigorous result on this 
model at positive temperature is probably in [22J and the best result up to date is in [9J; (see also Section l4~3l l. 

To describe some related results, we first introduce a different model, which, however, will turn out to be closely 
related to ours. 

Assume that the quantum particle interacts with random time-dependent impurities. That is, let V(x, t) be a 
real-valued random variable, for x € Z d , t £ R, with mean zero 

E [V(x, t)] = 0, (1.21) 

satisfying the Gaussian property 

E[V(x 2n ,t 2n )...V(x 1 ,t 1 )} II ®[V(x s ,t s )V(x r ,t r )} (1.22) 

pairings 7r (r,s)G7r 

E[V(x 2n+1 ,t 2n+1 )...V(x 1 ,t 1 )] = 0, (1.23) 

where a pairing it is a partition of {1, ... , 2n} into n pairs and the product is over these pairs (r, s). In addition, we 
assume that the correlation functions are invariant under translations in time and space, 

E [V(x, t)V(x', t')} = E [V{x -x',t- t')V(0, 0)] . (1.24) 

A time-dependent Hamiltonian is given by 

H x (t):=-A + X^2v{x,t)\x)(x\, (1.25) 
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and the dynamics U x is defined (almost surely) by 

^U t x = -\H x (t)U t \ U x = 1. (1.26) 

One can check that if we choose 

E [V(x, t)V(0, 0)] := S x> o f dg|^(g)| 3 (pi[ao(?)aS(g)]e- ifa '(9) +/ ^ [05(5)00(5)] e'Mrt), ( L2 7) 

(the RHS will be motivated in Section l2.3.2|l , then we have that 

E[U t x p(U t x Y] = Z x (p). (1.28) 

The reason for this equivalence is that both models share a "quasi-free", or, 'Gaussian property". (In the Hamil- 
tonian model, this is a consequence of the fact that the free reservoir Hamiltonian is quadratic in the creation and 
annihilation operators). Of course, it is not clear that the definition (|1.27) l makes sense. For example, the RHS could 
have an imaginary part, whereas the LHS is real. However, upon inspection of our proof, it becomes clear that 
whenever 

\E[V(x',t')V(x,t)}\ < S x , x ,ce- 9Rlt - t ' 1 , c<oo, 5R <oo, (1.29) 

then our proof (which assumes the same bound for the RHS of CGZj) carries over, and we can establish diffusive 
behavior for the averaged dynamics E [U x p(U x )*] . In fact, the locality in space (expressed by 8 XtX >) is not crucial, 
at all, but we do not pursue this generalization here. 

The case E [V(x' , t')V(x, t)] = 5 x . x 'S(t — t') has been treated in fl6l . While we were completing this paper, a 
preprint [15 [ appeared where diffusion is proven under the assumption that V(x,t) is an exponentially ergodic 
Markov process (not necessarily Gaussian) for each x. Preliminary results were obtained in [17 \ and [24J. One of 
the ultimate goals of these projects is to treat the case where V(x, t) = V(x) is time-independent and d = 3, i.e., 
E [V(x')y(x)] = S XjX /. This is the well-known Anderson model. 

Models in which the particle is coupled to a thermal reservoir are expected to be easier than the Anderson 
model, mainly because one expects that diffusion persists for large values of the coupling constant A, whereas the 
Anderson model has a phase transition, and the particle gets localized at large values of |A|. 

However, even for a particle coupled to a thermal reservoir in d = 3, our techniques fail, since this model 
would essentially correspond to one with E [V(x', t')V(x, t)] ~ \ x l x >\ x[\ x ~ x '\ > c \t ~ ( ror reservoir particles 
with dispersion relation uj(q) = c\q\). 

There are however results that establish diffusive behavior up to times of order \~( 2+s \ for some 5 > 0, even 
for the Anderson model, see fTTllTOl (a resulting lower bound for the localization length is proven in [5]). In fact, 
our technique employs results of the type proven in these references as an ingredient of the proof; see Section l43l 

We might add that we expect that the model treated in the present paper can also be analyzed using operator- 
theoretic techniques introduced for the study of return to equilibrium in open quantum systems, see e.g. fl4l |2) , 
and we are currently working on such a formulation. The technique used in the present paper is largely based on 

m 

1.4 Outline 

In Section |3 we introduce our model, making precise the description in the introduction. Then, in Section [3l we 
state our assumptions and main results with as few divagations as possible. Section [4] contains the main ideas of 
the paper and the plan of the proof. The technical parts of the proof are postponed to Section|5l which contains the 
proof of Theorem l4.41 and Section^ where one finds the proof of Theorem l4.5l . 

1.5 Acknowledgements 
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2 Model 



2.1 Conventions and notation 

Given a Hilbert space S ', we use the standard notation 

with (<?) = ^(<?) the bounded operators on and 



< 00 



1 < p < 00, 



HSIIp : = far [(5*^) p/2 j) 



i/p 



|5|| :=||5|| 



(2.1) 



(2.2) 



For bounded operators acting on S% v {<§), i.e. elements of @)(8§ v ($)), we use in general the calligraphic font: 
V, W,T, . . .. An operator A" e ^(<f ) determines an operator ad(X) € S{M V {S)) by 



adpQS := [X, 5] = - SX, S & M p (£). 
We will mainly use the case p = 2. The norm of operators in 38(382(&)) is defined by 

'|W(5)|| 2 



l|W|| 



sup 



\S\\ 



(2.3) 



(2.4) 



For vectors k e C d , we let !Rk, 3k denote the vectors (SRk 1 , . . . , 3?ft d ) and (3k 1 , . . . , 3K d ), respectively. The scalar 
product on C d is written as (•, ■) and the norm as \k\ :— (k, k). The scalar product on an infinite-dimensional 
Hilbert space S is written as (•,•), or, occasionally, as (•,)<?• All scalar products are defined to be linear in the 
second argument and anti-linear in the first one. 

We write T S (S') for the symmetric (bosonic) Fock space over the Hilbert space $ and we refer to |7] for defi- 
nitions and discussion. If uu is a self -adjoint operator on S , then its (self-adjoint) second quantization, dr s (w), is 
defined by 



dr s HSym( 



:= ^ Sym(0i ® . . . ® tufa ® 



0. 



(2.5) 



where Sym projects on the symmetric subspace and . 



2.2 The particle 

We set J^s = l 2 (Z d ) (the subscript S refers to 'system', as is customary in system-reservoir models). We define the 
one-dimensional projector l x on fflfe by 

(l x f){x) := 5 x , x ,f{x'), x, x' eZ d Je l 2 (Z d ). (2.6) 

We will often consider the space J#s i n its dual representation, i.e. as L 2 (T d ,dk) where T d is the d-dimensional 
torus, which is identified with L 2 ([— n, n] d ). We define the 'momentum' operator P as multiplication by k S T d , 
i.e., 

(P0)(fc) := k9(k), 6 E L 2 {T d ,dk). (2.7) 

Although P is well-defined as a bounded operator, it does not have nice properties; e.g., it is not true that [X 1 , P 3 ] = 
i5ij. Throughout the paper, we only use operators f(P) where / is periodic on M. d with period 2ir, i.e. a function 
on T d . We choose a periodic function e to be the dispersion law of the system. Although this is not essential, we 
require e to have inversion symmetry, i.e., 

e(k)=e(-k), keT 1 . (2.8) 

The Hamiltonian of our particle is given by 

H s :=s(P). (2.9) 
Our first assumption ensures that H$ is sufficiently regular. 
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Assumption 2.1 (Analyticity of system dynamics). The function e, defined originally on T d , extends to an analytic 
function in a strip of width S e > 0. That is, when viewed as aperiodic function on R d , e is analytic in (R+i[— S £ , S £ ]) d . 
Moreover, we assume that the function T d 3 k (y, Ve(/c)) does not vanish identically for any vector v 6 R d , v ^ 
0. 

The most natural choice for e satisfying Assumption 12.11 is e(k) — $2f =1 (2 ~ 2cos(fc : ')), which corresponds to 
— Hg being the discrete Laplacian. 

2.3 The reservoirs 

2.3.1 Reservoir spaces 

We consider an array of independent reservoirs. With each site x 6 Z d we associate a one-particle Hilbert space 
t) x (one can imagine that t) x = L 2 (R d )) with a positive one-particle Hamiltonian u x . The reservoir at x is now 
described by the Fock space r s (f) a; ) with Hamiltonian dr^u^). The full reservoir space is 

■= T s (® xezd t) x ) with Hamiltonian H K := ^ dT s (u x ). (2.10) 

We choose the different reservoir one-particle spaces to be isomorphic copies of a fixed space \) so that <p £ t) x 
is naturally identified with an element of \) x > that is also denoted by <p without further warning. Likewise, oj x is 
naturally identified with oj x i . Hence, if no confusion is possible we simply write f) and ui to denote the (one-particle) 
one-site space and the Hamiltonian, respectively. 

For ip 6 t), the operators a* x (<p)/a x (ip) stand for the creation / annihilation operators on the Fock space T s (t) x ) . By 
the embedding of t) x into (B ye z d ^yf these creation/ annihilation operators act on in a natural way. They satisfy 
the commutation relations 

[a x (v),a x ,(v')}=6 x<x ,{tp,<p , } l) , [a#(<p),c$(<pf)]=0 (2.11) 
where a* stands for either a* or a. 

2.3.2 Interaction and initial reservoir state 

We pick a 'structure factor ' <f> G f) and we choose the interaction between the system and the reservoir at site x to 
be given by 

l a where ¥ x (0) = a x {4>) + a* x {4>) (2.12) 

So far, we have not made any assumptions concerning u> and <f>, but their form will be restricted by Assumption 
12.21 in l|2.20t . The particle interacts with all reservoirs in a translation invariant way. Hence the total interaction 
Hamiltonian is given by 

#sr:= Yl on ^s®^R- (2-13) 

Next, we put the tools in place to describe the positive temperature reservoirs. Let ^ be the ^-algebra consisting 
of polynomials in a x (ip),a* x , (ip'), with tp, <p' £ f), x, x' S 2, d . We introduce the positive operator Xjg = (e 13 ^ — l) -1 on 
f); (3 should be thought of as the inverse temperature. 

We let be a quasi-free state defined on ^ . It is fully specified^ by 

1) Gauge-invariance 

^tefa)]=A&Mv)]=0. (2.14) 

2) Two-point correlation functions 

( pkK{v)<>{<p')] \ = s ( We?) o 

\ &\a*(<p)ax>W)] &[a x {y)a* x ,{ V >)] ) x ' x> \ {<p\(l + T fi W)) 

2 The reason why, in models like ours, it is enough to know the state on %' , has been explained in many places, e.g. 1 1 3 8 13 1. 
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3) Quasi-freeness, i.e. , the higher-point correlation functions are expressed in terms of the two-point function 
by 

/4[a#>i)...a#> 2 „)] = E II pi[4M)a*^ s )] (2.16) 

pairings tt (r,s)£7r 
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a*(^i)---af„+i(^2n+i) 



= 0. (2.17) 



where a pairing it is a partition of {1, ... , 2n} into n pairs and the product is over these pairs (r, s). 
A quantity that will play an important role in our analysis is the on-site-reservoir correlation function defined 

by 

= (0, l>e w <» + (0, (1 + Tp)e- H «4>), (2.18) 
It is useful to introduce i/j, the inverse Fourier transform of %l> 

4>(t) = -L / d£e*V(0- (2.19) 



2tt 

As is explained in Appendix A, tp is the (squared norm of) the effective structure factor . In particular, ?/>(£) > 0. 
The following assumption requires the reservoir to have exponential decay of correlations. 

Assumption 2.2. There is a decay rate gn > such that 



sup 

ten 



(|^(i)|e 9R| *') < oo. (2.20) 



We assume that $ ^ 0, or equivalently ?A ^ 0. 



The assumption that =k ensures that the particle interacts effectively with the fields describing the reservoirs. 
In Appendix A, we discuss examples of reservoirs that satisfy Assumption ^. 21 provided that (3 < oo. 

2.4 The dynamics 

Consider the zero-temperature Hilbert space J?s ® J^k- The Hamiltonian (with coupling constant A) is formally 
defined by 

H x := H S +H R + XH SR . (2.21) 

This operator generates the zero-temperature dynamics. However, we need to consider the dynamics at positive 
temperature. In particular, we must understand the reduced positive-temperature dynamics of the system S after 
the reservoir degrees of freedom have been traced out. 

By a slight abuse of notation, we use p R to denote the conditional expectation from £§(J$?s ® to 33(J%%) given 

by 

pi(S® R) := Sp R (R), Se,^(,mi),R&^ (2.22) 

where p R {R) is defined through ( |2.15i2.17> . 

Formally, the reduced dynamics in the Heisenberg picture is given by 

Z^* (S) := p R [c itH > (S (g) 1) e- itH "] (2.23) 

whenever the RHS is well-defined. 

A mathematically precise definition of the reduced dynamics is the subject of the next lemma. 
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Lemma 2.3. Suppose that Assumption \2.2\ (see M.20\ ) holds and define 

H SR (t) := ® *x(e i4 ^) iwtfi l a (f) :=e ims l a e- itHs . (2.24) 

T7ze serie^ 

Z X '*(S) := ^ (iA)" / dti . ..dt n /4 [ad(JTsn(ti)) • . . ad(tf S R(U)) e itad(ffs) (S ® 1) ] (2.25) 



0<ti<...<t„<t 



zs well-defined for any A,t e R awd arbitrary S £ i.e., the RHS converges absolutely in the norm of and 

Z x '* has the expected properties, namely 

Z t A '*(l) = l, \\Z^(S)\\<\\S\\. (2.26) 

One can prove this lemma (under less restrictive conditions than those in Assumption l2.2b by direct estimates of 
the RHS of 02.251 1. For this purpose, the estimates given in the present paper amply suffice. However, one can also 
define the system-reservoir dynamics as a dynamical system on a Von Neumann algebra through the Araki- Woods 
representation. This is the usual approach in the mathematical physics literature; see e.g. I8 H141IT31 . 

Finally, we define Z t A : 8§i{J#s) — > the reduced dynamics in the Schrodinger picture, by duality, i.e., 

Tr[ Ps Z^*(S)} = Tr[Z x (p s )S], S G ^(Jf S ),p S € #i (.*§). (2.27) 

We could also have started by defining the full initial state psn of the total system consisting of the particle and 
reservoirs as the positive, normalized functional 

PSR~PS®P^ on 39(3%)®*, (2.28) 

where we abuse notation by employing the same symbol p$ for both the density operator (a positive element of 
^i(J0§)) and the state it determines on SS(J^ S ), i.e., 

p s [S] :=Tr[p s S], Se» (2.29) 

Then, 

PSR [e itH > (S S 1) e-"^] = Tr[Z t \p s )S]. (2.30) 

In what follows, we simply write p for pg. For convenience, we treat p as an element of the Hilbert space ^(^s)/ 
which is justified since ^i(^) C 3§ 2 {<%k)- 

3 Result 

We now state our main results. Recall that the position operator X on l 2 (Z d ) is given by 

(Xf)(x) = xf(x), xeZ d ,fel 2 (Z d ). (3.1) 

For k € C d , we define 

J K S :^e-^ x) Se-^ x \ SeJ(Jf s ). (3.2) 
Note that J K is unbounded if k £ R d . We choose an initial state p e 8$i(J%s) satisfying 

p > 0, Tr[p] = 1 \\J K p\\ 2 < oo, (3.3) 

for k in some open neighborhood of G C d . 

Our first result says that the momentum distribution of the particle tends to a stationary distribution exponen- 
tially fast. 



3 In fact, one needs to do things more carefully, since -HsrW t ■ A possible solution is to define the cut-off interaction -Hs-r,a(*) = 
J2 X £A ® ^(e 1 '" 1 / 1 ), for some finite subset A C Z d , and to show that one can take the limit A / Z, d in the expression analogous to <2.25> . 
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Theorem 3.1. [Equipartition Theorem] Suppose that Assumption \2.1\ (see Section \2.2\ and Assumption \2.2\ (see i2.20i ) 
hold, and let p satisfy condition ||3.3|I . There are positive constants \ > and g > such that for < |A| < Ao, there is a 
function C° € L 2 (T d ) satisfying 

Tr[8(P)Z>>(p)} = (0,t x ) L 2 m +O(\\0\\ 2 e- x29t ), ast / oo, for any 8 = 6 e L°°(T d ), (3.4) 

and 

-/3e(fe) 

m = j^ m + 0(X% AXO. (3.5) 

The decay rate X 2 g is strictly smaller than gn, introduced in (|2.20l l. 

Define a probability density p x depending on the initial state p G SS\ (J%s ) by 

p x (x) :=Tv[l x Z x {p)] . (3.6) 

It is easy to see that 

p x (x) > 0, ^{x) = Tr[p] = 1. (3.7) 

We claim that the particle exhibits a diffusive motion. This is the content of the next result. 

Theorem 3.2. [Diffusion] Under the same assumptions as in Theorem \3~T\ the following holds. Let the initial state p satisfy 
condition l|3.3b and let p x be as defined in (|3.6l l. There is a positive constant Ao such that, for < |A| < Ao, 

Y, /i* A (z)e~^ (9 ' x) e -^ D >o\ q e K d (3.8) 

ztfere f/ze diffusion matrix D x is positive-definite (i.e., has strictly positive eigenvalues), and 

D\ = A" 2 (£> kin + 0(A 2 )) , A \ 0, (3.9) 

with £>kin a \-independent positive-definite matrix introduced in Section \4~3\ 

We refer to Section [TJ for an explanation of the connection between this result and diffusion in the physicists' 
sense. We close this section with some remarks concerning possible extensions of our results. 

Remark 3.3. Our proof of Theorem \3\2\ actually gives a stronger result. Assume the nth moments of the initial distribution 
are bounded, or, equivalently, 

<jn 2_, l 1 o{ x ) e ~ lt ' q,x " > is n times differentiable. (3.10) 

Then the rescaled nth moments converge to the nth moments of the limiting distribution, or equivalently, the derivatives of 
nth order of 

rf(x>~^ iq ' X) (3-11) 

converge, ast /* oo, to the derivatives ofe~( q ' Dxq \ For n = 2, this implies ((1.191 . Note that the condition 03.101 1 is a weaker 
assumption than ((3.3b ; in fact, (|3.3b implies that d3-101 > is a real-analytic function. 

Remark 3.4. By the same technique as employed in our proofs, one can show that correlations decay rapidly in time. As 
explained in the introduction, this rapid decay provides an intuitive explanation why the particle motion is diffusive. 
Define the particle velocity operator by 

V{t):=ie ltHx [H x ,X]e- itH ^ (3.12) 

and observe that 

V(0) = i[H x ,X] = i[H s ,X] = (Ve)(P). (3.13) 

Suppose that Assumptions \2.1\ and \2.2\ hold and let p — ps satisfy condition i3.3i . 

By reasoning similar to that in Lemma \2~3[ one can define the velocity-velocity correlation function psr [V (t\)V (£2)]- Let 
the coupling strength A and the positive constant g be as in Theorem \3~T\ Then, for all < t\, t 2 < oo, 

|psR[V(ti)V(t 2 )] I < ce- x2alt2 - tl1 , forsomec<oo. (3.14) 
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Remark 3.5. The condition that the particle dispersion satisfies e(k) = e(— k) is not really necessary for our results to hold. 
If one did not impose this condition, the particle could have a drift velocity v^t given by 

«dr:=(Ve,CS>, (3.15) 

and the particle motion would still be diffusive, but one would now consider the "random variable" -^(xt — VdA), instead of 
-^x t . In other words, in (|3.8b , one would have to replace 

J2 ^ t {x)e-irM^ by rf(x>-^ {9 ' (x - Vdrt)) - (3-16) 

Similarly, in eq. l|3.14l l, one would have to replace V(t) by V(t) — v<i T , 

4 Discussion and outline of the proof 
4.1 Translation invariance 

Consider the space of Hilbert-Schmidt operators S§ 2 {^s) ~ & 2 (l 2 {Z d )) ~ L 2 {T d x T d , dfcidfc 2 ), and define 

S(h,k 2 ) := £ S{x 1 ,x 2 )e-^ k ^^\ 5€J 2 (i 2 (Z d )). (4.1) 

In what follows, we simply write S for S. To deal conveniently with the translation invariance in our model, we 
make the change of variables 

fci + fc 2 , , „. 

k= — - — , p = k 1 -k 2 , (4.2) 

and, for a.e. p S T d , we obtain a well-defined function S p G L 2 (T d ) by putting 

(S p )(fc):=S(fc + f,fc-f). (4.3) 

This follows from the fact that the Hilbert space SS 2 {Mk) ~ L 2 (T d x T d , dkidk 2 ) can be represented as a direct 
integral 

= / d P Jf p , S= I d P S p , (4.4) 
where each 'fiber space' .yf v is naturally identified with L 2 (T d ). Let T z ,z G Z d , be the lattice translation 

(T z S)(xi,x 2 ) := S(xi + z,x 2 + z), S € &(J%), (4.5) 

or, equivalently, 

(T z S) p (k)=e i(p ^S p , Sef(Jf s ). (4.6) 
Since Hx and p# are translation invariant, it follows that 

T_ z Z t A T 2 = Z 4 A . (4.7) 

Let W € &(& 2 (J?s)) be translation invariant in the sense of eq. j47fr . i.e., T_ 2 WT 2 = W. Then it follows that, in 
the representation defined by (|4.4) l. W acts diagonally inp, i.e. (>VS) P depends only on S p , and we define W p by 

(WS) P = W P S P . (4.8) 

For the sake of clarity, we give an explicit expression for W p . Define the kernel W(x, y; x' , y') by 

(WS)(x',y')= W(x,y,x',y')S(x,y), x',y'eZ d . (4.9) 

x,y& d 
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Translation invariance is expressed by 

W(x,y;x',y') = W(x + z,y + z:x + z,y' + z), z£Z d , (4.10) 
and, as an integral kernel, W p £ S3(L 2 (T d )) is given by 

W P (k',k)= Yl jt k *-v)-WJ-v%iW x W-(*^»W(x,y;x!,i/). (4.11) 

x,y,x',y' £ Z d 
x + y + x' + y' = 

Next, we state an easy lemma. 
Lemma 4.1. Let S £ SS\ (<3#s ) . Then, S p , as defined in i43\ , is well-defined as a function in L 1 (T d ) for every p, and 

Tr[J p S] = e- ipx S(x,x) = (1,S P ). (4.12) 

x£Z d 

where 1 £ L 2 (T d ) n L°°(T d ) is the constant function with value l(fc) = 1. Assume, moreover, that there is a constant 5 > 
such that 

\\J K S\\ 2 < oo for |9f/s| < S, (4.13) 
then the function p i— > S p £ L 2 (T d ) has a bounded-analytic extension to the strip |SJp| < S. 

The first statement of the lemma follows from the singular-value decomposition for trace-class operators and 
standard properties of the Fourier transform. In fact, the correct statement asserts that one can choose S p such that 
04.121 1 holds. Indeed, one can change the value of the kernel S(ki , k 2 ) on the line k\ — &2 = P without changing the 
operator S, and hence S p in (|4.12) l can not be defined via l|4.3|l for all p, but only for almost all p. 

The second statement of Lemma 14.11 is the well-known relation between exponential decay of functions and 
analyticity of their Fourier transforms. Since we will always demand the initial density matrix po to be such that 
1 1 J7k Po 1 1 2 is finite for k in a complex domain, we will mainly need the second statement of Lemma |4~T1 

4.2 Return to equilibrium inside the fibers 

The main idea of our proof is that the reduced evolution in the 'low momentum fibers', (2^) p , for p near 0, has an 
invariant state to which every well-localized initial state relaxes exponentially fast. 

Recalling that Hs = e{P) and that the system is weakly coupled to a heat bath at inverse temperature /3, we 
expect that, in an appropriate sense, and for arbitrary initial states p £ £$i(Hs), 

tjoo £\P) 

We observe that or^ ei - p ^ ^ SS\{H^), hence (|4.141 cannot hold in norm (in other words, Z(/3) = oo). One way to 
interpret (14.141 is that it gives the correct asymptotic expectation value of functions of the momentum, and that is 
exactly what Theorem 13.11 states. 

For every p satisfying | |3.3l l, we have that 

Tv[d(P)Z t x (p)} = (9, (Z t x p) ), 9 £ L°°(T d ), (4.15) 

by applying Lemma [4 . 1 1 with S := 9(P)Z x (p). Hence, we should apparently attempt to prove 'return to equilib- 
rium' for the evolution (Z x ) on L 2 (T d ). 

The dynamics in the fibers corresponding to small values of p provides information on the diffusive character 
of the system. The probability density p x (x) corresponding to some initial state p is defined as in l|3.6t . By Lemma 

^rt(x)e- i M=J2( z tP)fa*)e~ i(p '* ) = f /k{Z x p){k + V -,k- P -) = {l,{Z x p) p ). (4.16) 

x£Z d x& d Jd 
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To establish diffusion, it suffices to show that, for A fixed and for p in a neighborhood of G T d , 

<l,(Z ( V)p) =c t( --^ p ' D ^ + °( p2 »(l + o(t°) + o( P )), t/*oo,p\0. (4.17) 

for some positive-definite matrix D\. Indeed, by | |4.16| |, Theorem |3.2| follows from (|4.17l l by taking p = Thus, in 

order to prove Theorem |3.21 we are led to study the long-time asymptotics of the evolution (Z^) p , for small p. 

However, as our approach is perturbative in A, expression l|4.17l l is not a good starting point, since (p, D\p) = 
0(A~ 2 ), for fixed p (as can be seen from the statement of Theorem 13.2b , and hence one cannot perturb around 
(p, D\p)\ x _ Q . The way out of this difficulty is to set up the perturbation on a scale where the diffusion constant is 
finite (this will turn out to be the kinetic scale), or, in other words, to take the p-neighborhood in l|4.17|l to shrink, as 
A \ 0. Since A approaches 0, one must wait a time of order A -2 , before one sees the effect of the interaction. Since, 
between collisions, the velocity of the free particle is unaffected, it travels a distance of order A~ 2 . This means that 
when both space and time are measured in units of A~ 2 ; 

x = \- 2 i x , t = \~ 2 tx, (4.18) 

we expect a diffusion constant D\ ~ {x ~ x> of order 0(1). This is consistent with the fact that D\ ~ ^ is of order 
A~ 2 . The limit D\\o is the diffusion constant in the kinetic limit, as outlined in the next section. 

4.3 The kinetic limit 

To control the asymptotics of the effective time-evolution (Z^) p , we compare it with the corresponding evolution 
in the kinetic limit, which is the limit approached when microscopic space and time are taken to be A~ 2 x, A~ 2 <, 
respectively, and the coupling strength A — * 0; as announced in the previous section. It has been proven in [9 [ 
(for models with only one thermal reservoir) that, in this limit, the dynamics is described by a linear Boltzmann 
equation. 

Our variant of this result is described below. 
4.3.1 Convergence to a linear Boltzmann equation 

The effective reservoir structure factor ijj has been defined in l l2.18H2.19l l. For convenience, we introduce a positive 
function r(-, •), with 

r(k, k') := ip[e(k') - e(k)] > 0. (4.19) 
For k £ M. d , we define a bounded linear operator, M K , on L 2 (T d ) by 

(M K e)(k) :=i(K,Ve)(A:)0(fc) + f dk' [r(k' , k)6(k') - r(k, k')0(k)] , 9 e L 2 (T d ), (4.20) 

where (k, Ve)(fc) stands for the scalar product in C d of k and We(k). The operator M K has a straightforward 
interpretation: Consider a classical particle whose states are specified by a position x E M. d and a 'momentum' 
k € T d . The momentum k evolves according to a Poisson process with a rate r(k, k') for the transition from state k 
to k'. Between two momentum jumps, the particle moves freely, with speed given by (Ve)(/c). The translation of 
this picture into a mathematical statement is as follows: The state-space distribution of the classical particle at time 
t is given by a probability density v t {-, ■ ) on R d x T d ; (v(x, k) > and J dxdk v t (x, k) = 1). Then 

^-u t (x ) k) = (V k e,V x p t )(x,k)+ f dk' [r{k\k)y t {x,k')-r{k,k')v t {x,k)]- (4-21) 
at j T d 

One checks that 

v?{k) := (27r)- d / 2 f dx e-' l( - K ' x K t (x, k) (4.22) 

JR d 

satisfies an evolution equation generated by M K ; 

= M K i>t . (4.23) 
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We claim that the rates r(k,k') satisfy the identity 

r(k, k') = r(k', k)e~^ E(k ">- e( - k ^ , (4.24) 

known as the detailed balance condition in the context of Markov processes. It is a direct consequence of the KMS- 
condition for the reservoirs. In our context, it is easily derived from {2.15) . The detailed balance condition implies 
that 

e -(3e(k) 

M°CL = 0, where CL(*0 = (4-25) 

In the language of Markov processes, £ki n is a stationary state. 

The relevance of M K is that it describes the evolution Z^_ 2t in the fiber indexed by X 2 k in the limit A \ 0. 
Moreover, the convergence of the fiber dynamics (Z*_ 2t )\2 K holds even after analytic continuation to complex k. 
One can prove the following result 

Proposition 4.2. Assume Assumptions \2.1\ and \2.2\ Then, for \Qn\ sufficiently small and < t < oo, 



c 



LM" 



— ► 0. (4.26) 

A\,0 



where the norm is the operator norm on L 2 (T d ). 



We do not prove this Proposition (which is not needed for the proof of our results). In fact, the proof is based on 
the same reasoning as in Section[6] Of course, one can also express Proposition 22] in terms of the rescaled Wigner 
function, as is done in EKE]- Indeed, setting 

a?(k) := lim (Z^p) (k + X 2 ^k- A 2 ^) = lim AZ^ H p) x , K {k), (4.27) 

one obtains from Proposition l4.2l that a" (k) satisfies the evolution equation (|4.23l l. (It would thus be justified to call 
a.t(k) simply 0^(k)). Its inverse Fourier transform 

a t (x,k) = (2Tr)- d / 2 { dKe iiK ' x) a?(k) (4.28) 



is a probability density on K d x T d and satisfies (|4.21) l with initial condition cto(x, k) = S(x)p(k, k). 

We state another useful consequence of Proposition 14.21 Recall that the probability density pt(-) has been de- 
fined in (|3.6l l, for any initial state p. Taking the scalar product with 1 S L 2 (T d ) on both sides of l!4.27t and using 
II4.16I I, we obtain that 

J2 e-** M A- H (x) dkaUk). (4.29) 

As outlined in Section l4~2l the t f oo asymptotics of the LHS of (|4.291 contains information on the diffusive 
behavior of the particle. In the next section we discuss the t f oo asymptotics of the RHS of l!4.29l l. 



4.3.2 Diffusive behavior of solutions of the Boltzmann equation 

To realize that the Boltzmann equation describes diffusion, one studies the spectral properties of M K , for small 
k. We state a crucial result, Theorem 14.31 and we refer the reader to |6| for complete proofs and a more extended 
discussion of quantum dissipative evolutions. 

Theorem 4.3. Suppose that Assumptions\2A\and\2j\hold, and let M K e £g(L 2 (Y d )) be defined as in (4.20)1 . 

Then there is a positive constant 5un such that the operator M K , with \k\ < S^, has a simple eigenvalue, /ki n («0/ 
separated from the rest of the spectrum by a gap, 

dist(5R/ kin ( K ), m) =: g k in > (4.30) 
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where 

n := U w<4kin (s P M K \ {/ k in(K)}) , anrf 3?fi < 0. (4.31) 
The eigenvalue / kin and its associated eigenvector Q in and spectral projection P£ in are analytic in k and 

/kin(«) =-<l,(«, V £ )(M )" 1 ( K ,Ve)C k ° in ) + O( K 3 ), «\0, (4.32) 
iw/zere Ve denotes the operator that acts by multiplication with the function Ve on T d . T/ze diffusion matrix, D^, defined by 

(Adn)^' := -^jAin(«)L =0 , i,i = 1, •• -d, (4.33) 
ftfls real entries and is positive-definite. 
Sketch of proof. 

We write M° =K + T with 



(Kff){k)= I dk'r(k',k)6(k'), (T9)[k) = - [ dk'r(k,k') ) d(k), 9 e L (T , dfc). (4.34) 
Notice that T is a multiplication operator with spectrum 



spT=|-y dk'r[k,k') \ k e T d | , r(fc,fc')>0. (4.35) 

The operators K and T are sometimes referred to as the gain and loss terms in the Boltzmann equation. Assump- 
tions |2J] and [272] imply that the functions ip and e are real-analytic in k, and hence r(-, ■) is real-analytic in both 
variables. It follows that K is a compact operator on L 2 (T d ) and, since we assumed that tp ^ 0, we have that 
sup^spT < 0. By Weyl's theorem on the stability of the essential spectrum (see e.g. p. 101 of [18]), we deduce 
that the spectrum of M° in the region > sup SftspT consists of isolated eigenvalues of finite multiplicity. From 
the pointwise positivity of r(-, •), the Perron-Frobenius theorem and from the fact that M° generates a contractive 
semigroup on L 1 (T d ) we then conclude that the eigenvalue of M° is simple and that it is separated by a gap from 
the rest of the spectrum. The spectral projection P kin is explicitly given by 

P k V = (MRkV 9eL 2 (T d ) (4.36) 

with as in J4.25|l . The analyticity of /kin(t) and ££ in is proven with the help of analytic perturbation theory. 
Using the assumption that e(k) = e(~k), we check that 

PLVePL = 0- (4.37) 
Employing explicit expressions of second order perturbation theory, we obtain formula J4.32|) as a consequence of 
the fact t hat M K - M° = Ve) and @37|> - 

Since /ki n (/t) = /kin(— k)/ it follows that the matrix -Dkin has real entries. The positive-definiteness of £>kin is 
established as follows. Consider the bounded operator 

(W9)(k) = e^ e(fc) 6>(fc), 9 e L 2 (T d ), (4.38) 

and notice that M := W~ 1 M°W is a self-adjoint operator on L 2 (T d ), in particular £ := WCkin = W~ x l, (i.e., the 
left and right eigenvector corresponding to the eigenvalue are identical). For n e M. d , we can rewrite l|4.33l l as 

(k, Ad^) = <(«,Ve)C, Ve)C>. (4.39) 

By Assumption 12.11 the function fc i— > (k, Ve(/c)) does not vanish identically on T d (for k ^= 0). Hence, by the 
spectral theorem applied to M, expression J4.39& is strictly positive. □ 

Let Vt{k) be a solution of the evolution equation | |4.23|| for k in some neighborhood of in C d . Using Theorem 
I4.3l and reasoning similar to that in Section l4~2l it follows that 

dku?(k) — ► e -^ (9 ' Dk '» 9) , q eR d . (4.40) 

Hence a solution v t {x, k) of the Boltzmann equation (|4.21l l behaves diffusively, with diffusion tensor -Dkin. 
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4.4 Perturbation around the kinetic limit 



Up to now, we have seen that, in the kinetic limit, the particle motion is described by a linear Boltzmann equa- 
tion. Since solutions of the linear Boltzmann equation behave diffusively for large times (as is essentially stated in 
Theorem |4.3l l, we can associate a diffusion constant to our model. Indeed, by d4.29l l and ( 14.401 ), 

lim lim ^x-H( x ) e ^ (q ' X) =e-5fe' D ^9) (4.41) 

However, (14.411 does not give information on the long-time asymptotics of our system for small, but fixed |A| > 0. 
The least one would wish for is to be able to exchange the order of limits in d4-41b , and, indeed, Theorem |3.2| states 
that one can do so without affecting the RHS. We stress this point, because it is an improvement of our paper when 
compared to most earlier results on diffusion. 

Since we have learned that (Z^_ 2t )\2 K has a well-defined limit, e tM , as A \ 0, (see Proposition l4.2|l , it is natural 
to expand (Z^_ 2( )a 2 k around this limit, in such a way that we can take t f oo. We perform the expansion on the 
Laplace transform of Z t A , 

Ux(z) := / dte~ tz Z t x . (4.42) 



Theorem l4.4l below summarizes the result of our expansion. Loosely speaking, a key consequence of this theorem 
is the fact that, in the fibers indexed by A 2 k, one has that 

{n x (z)) X 2 K = (z- \ 2 M K - A(z, A, k))~\ (4.43) 

where the operator A(z, A, k) is "small" compared to \ 2 M K . 

Theorem 4.4. Suppose that Assumptions \2.1\ and 12.21 in Section^ hold. Then, there are operators L(z) and TZ° x x (z) in 
£M{&2{3%&)) such that the following statements hold. 

1) For (z,A)eCxl satisfying 3?z > ||A 2 £(z) + K e f(z)\\, 

TZ x (z) = (z- ad(Lffs) - A 2 £(z) - nf(z))-\ (4.44) 

2) The operators C(z) and lZ e ^(z) have the following properties: There are positive constants 8[,5 2 ,g' > such that 

J K C{z)J- K , J K 1lT(z)J_ K (4.45) 
are analytic in the variables (z, n) e C x C d in the region defined by \n\ < S[, $lz > —g',\X\ < 5' 2 . Moreover, 

sup \\J K C(z)J- K \\ = 0(1), A\0 (4.46) 

|k|<5J,5Rz>-9' 

sup \\J K Kf(z)J- K \\ -0(A 4 ), A\0, (4.47) 

|k|<<5J,sRz>-9' 

where \\ ■ \\ refers to the operator norm on ^(^(^s)) (us in (I2.4I >). 

3) Let M K be defined as in Section \4~3\ Then 

\\(nd(iH s ) + \ 2 £(0)) x2K -\ 2 M K \\=O(\ 4 K 2 ) + O{\ 4 K), A 2 k\0,A\0. (4.48) 

The proof of Theorem 14.41 is the subject of Section [5] From that proof, it becomes clear that g' can be chosen to 
be any fraction of gR by making 5 [ and 5' 2 small enough. 

From Theorem l4.4[ one obtains our main result by using Theorem l4.3l and standard analytic perturbation theory. 
More precisely, we prove the following theorem. 
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Theorem 4.5. Suppose that Assumptions \2.1\ and \2.2\ in Section\2\hold. Then, there are positive constants 61,82, g > such 
that, for (A, n) e K x C d and \k\ < 8\, < | A| < 82, there is a rank 1 operator P Xk and a function /(A, n) satisfying 

\\(Z t x ) X 2 K -e tf( - x - K '>P x - R \\^0(e t( - f ^- x2 ^), t/00 (4.49) 

and 

\\P x > K -PU\=0{\ 2 ), |/(A, K )-A 2 /kinWI -0(A 4 ), A\0 (4.50) 
Moreover, P X > K and /(A, n) are analytic in n g C d in the region defined by \n\ < Si, |A| < S 2 - 

By making 82 small enough, the constant g can be chosen to be any fraction of gkin and 81 can be chosen to be 
given by Ski n , with g^ n , <5kin as in Theorem |4.3t . 

Theorems 13.11 (Equipartition Theorem) and 13.21 (Diffusion ) then follow as discussed in Section l4~2l We briefly 
recapitulate our reasoning. 



Proof of Theorems \3.1\ and 13.21 

We first prove TheoremEO Using (|415l and Theorem|431 we write, for 9 = 6 e L°°(T d ), 

Tr[6(P)Z x p] = (6, (Z x p)o) = (9, e tf( - X '^P X '°p ) + 0( e *C/^°)- A2 s)). (4.51) 
Since Z x p has trace 1 (it is a density matrix) for all t > 0, we deduce that /(A, 0) = and, setting 9=1, 

(1,P A '%) =1. (4.52) 
The fact that P x,a is a rank 1 operator (by Theorem |4.5l l implies, together with | |4.52| |, that, 

P x >° v = (1(l,ri), foranyr/e i 2 (T d ), (4.53) 
for some € L 2 (T d ) which satisfies (1, C") = 1- Theorem 13.11 follows. 

We define the diffusion matrix by 

{Dx)l ' 1:= - X ^d^ f{X ' K) ^=o' '-' '••••' / - (4 ' 54) 



From | |4.12|| , with S := Z x p, we conclude that /(A, k) = /(A, —re), and hence the matrix D\ has real entries. Positive- 
definiteness of D\ follows then from positive-definiteness of -Dkin/ for A small enough. Using Theorem l4.51 we find 
that 

J2 e-^ x) p x (x) = (l,(Z x p)^) (4.55) 

= {1,{Z X P )^ K ), with K = A- 2 ^=, q eR d 

= (l,e t ^ x ^P x ^ PX 2 K )(l + 0(c-a t )), asi/oo 

= (l,e~ t(x4 ^ K ' D ^ )+o ^P x ' p Q )(l + O(K))(l + O(e- 9t )), a SK \0 

= (l,e-^' l ' D ^ +o ^P x > pQ)(l + O(K)){l+O(e- 9t )), 

which proves Theorem l3.2l upon using (1, P x >°p ) = 1 and n = A~ 2 ^=. □ 
Remark l3~3l follows by standard reasoning, using the following facts 

1) The family of operators 

{Z X ) X 2 K - e tf(*,K)p*,K (4.5 6 ) 
is analytic in k in a neighborhood of € C d and bounded by a constant independent of n and t. 

2) The function /(A, k) and the rank 1 operator P X ' K are analytic in k in a neighborhood of € C d 

This is related to the general fact that the central limit theorem follows from the existence and analyticity of the 
large deviation generating function, as described in [4J. Indeed, k /(A, n) can be viewed as the large deviation 
generating function corresponding to the family of random variables x t , t > 0, as defined in 111.171 1. 
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5 Dyson expansion and proof of Theorem 14^4 

To construct a "polymer model", we first write a Dyson expansion for Z t \ 



5.1 Dyson Expansion 

In this section, we set up a convenient notation to handle the Dyson expansion, which has been introduced in 
Lemma |Z3l Define the unitary group lit on ^ 2 (<^t) by 

U t S := e~ itHs Se itHs , S6f 2 (4), (5.1) 

and the operators I Xt i, with x E 7L d and I G {L, R} (L, R stand for "left" and "right"), as 

f i 1,5 if I = L 
I x iS := I , (5.2) 

1 -i Sl x if l=R. y ' 

Let P„ be the set of partitions it of the integers 1, . . . , 2n into n pairs. We write (r, s) E it if (r, s) is one of 
these pairs, with the convention that r < s. Note that the same notation was already used in i ll .221 1 and in (|2.16l l. 
Elements inR 2 ™, (Z d ) 2n , {L, R} 2n are denoted by t, x,l, with U,xi, k their respective components, for i = 1, . . . , 2n. 
We evaluate d2T25t by using j2/15b and j27T6] >- (f27T7b : 

(II E E Ut>&QUt-t 2n Z X2n ,i 2n ■ ■ .l X2 , h U t2 - tl l Xuh U tl (5.3) 



0<tl...<ta„<i 1=1 £i " G ' P " 



where 



sun TT x J ^ s ~ tr ) ^ = L ' « /n 



(r,s)^7r 

and, for n = 0, the integral in j5.3|l is meant to be equal to U t . 

We introduce some more terminology, extending the above definition of pairings. It will be helpful in classifying 
the pairings. 

Definition 5.1. 1) Let S n be the set of sets of n pairs of (distinct) natural numbers. More concretely, for each 
a E £ n , we can write 

a = {(ri,si),...,(r n ,s n )} ,n,Si E N, (5.5) 

for natural numbers r,-, s,-, i — 1, . . . , n which are all distinct. By convention, rt < Si,i — 1, . . . ,n and r$ < 
r i+ i,i = 1, . . . , n — 1. If ci G S ni and cr 2 G S„ 2 , we write o\ < a 2 whenever all elements of the pairs (rj, sj) 
in a\ are smaller than all elements of the pairs (r 2 ,s 2 ) in a 2l i.e., 

s-<r 2 , i = 1, . . . ,m, j = 1, . . . ,?i 2 . (5.6) 



2) Recall the definition of 7-Vi, the set of pairings with n pairs. Obviously V n C S n , and a E X„ belongs to P n 
whenever U" =1 {rj, s^} = {1, . . . , 2n} . Further, with any a E S„, we associate the unique pairing 7r E V n for 
which there is a monotone increasing function q on {1, ... , 2n} such that 

(M) E tt^ (q(i),q(j)) E a. (5.7) 

3) We set := ^ n >\P n and write |-7r] = n whenever tt e V n . 

4) We call a E S n irreducible (Notation: irr. ) whenever there are no two sets o\ E E ni , 02 € S„ 2 , ni + n 2 = n 
such that (7 = (7 1 U (72 and ci < 02- For any cr G E n that is not irreducible, we can thus find partitioning 
subsets cti, . . . ,a m (U™ x Oi = a) such that cr^i,. .., m are irreducible and <7j < er s ; + i for i = 1, . . . ,m — 1. 
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Figure 5.1: Graphical representation of a term contributing to the RHS of |53) with tt = {(1, 3), (2, 4), (5, 8), (6, 10), (7, 11), (9, 12)} G V&. 
The times t, correspond to the position of the points on the horizontal axis. 

Starting from this graphical representation, we can reconstruct the corresponding term in <5.3t - an operator on 3&2(J%))- as follows 

• To each straight line between the points (tj, Xi,li) and (tj+i, Xj+i, h+i), one associates the operators t4 i+1 — v 

• To each point (t;, Xi, U), one associates the operator X 2 I x . t i., defined in 15.21 . 

• To each curved line between the points (t r , x r ,l r ) and (t a , x a , l s ), with r < s, we associate the factor 



Rules like these are commonly called "Feynman rules" by physicists. 

5) Consider some it E V and its partitioning into irreducible subsets <j\, . . . , a m , as defined above. By \5.7\ , we 
can associate to each of the Oi a unique 7Tj in V. We call the set (ni, . . . ,n m ) of pairings, obtained in this way 
the decomposition of tt into irreducible components. 

6) For each n 6 N, we define a distinguished pairing tt e V n , which is called the minimally irreducible pairing 
(Notation: min.irr. ). For n > 2, this minimally irreducible pairing is given by 



For n = 1 and n = 2, the minimally irreducible pairing is defined to be (1, 2) and {(1, 3), (2, 4)} respectively. 
Intuitively, the minimally irreducible pairing in "P n is characterized by the fact that if one removes any pair, 
other than the pair with r = 1 or s = 2n, the resulting pairing is no longer irreducible. 




(ri,si) = (1,3), (r„,s. 



n 



) = (2n-2,2n), {r i+1 , s i+1 ) = (2i, 2i + 3), for i=l,...,n-2. 



(5.8) 



For an irreducible pairing tt E V n , we introduce (using the same conventions as in 115.31 1, ||5.4| |), 



2n-l 




t-i 2 




(5.9) 



0=tl<...<*2,i=t 



i=2 x,l 
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Figure 5.2: Graphical representation of a pairing it £ Pg. The pair (r, s) belongs to 7r whenever the natural numbers r, s are connected by an 
arc. This type of diagrams differs from those of Figure l5Tl in that we don't keep track of the tj -coordinates, but only of the topological structure 
of the pairings. Below is the decomposition of n into irreducible components. 




TTl 7T2 7T3 7T4 



Figure 5.3: The irreducible components wi, ^2,773,774. Explicitly, 7ri = 7r 2 = {(1,2)}, 7r 3 = {(1, 6), (2, 3), (4, 7), (5, 8)} and 7r 4 = 
{(1, 3), (2, 5), (4, 6)}. The pairings tti, 7T2 and 7T4 are minimally irreducible, whereas 773 is not. Indeed, one can remove the pair (4, 7) from 7T3 
without destroying the irreducibility. 

We can now rewrite d5.3b as a sum over collections of irreducible pairings; 



dti . . . dt 2r , 



m 6 Z+ 0<t!...<t 2m <t 



^^fe.Vfe.-^W . . .Wta-taVfa-t^TriJWt,. (5.10) 



7Ti , . . . , 7T Tn £ T 3 

7Ti , . . . , 7r m irr. 



To obtain this last expression, we decompose each pairing tt in (|5.3b into its irreducible components tt\, . . . , ir m , and 
we made use of a simple factorization property of l|5.3b . The term on the RHS of l!5.10b corresponding to m = is 
understood to be equal to lit ■ 

In expression (|5.10K we view the pairings tt^ with > 2 as excitations. If ]7r»| = 1, for alH = 1, . . . ,m, the 
corresponding term in (|5.10b is called a ladder diagram. These ladder diagrams provide the leading contribution 
to the dynamics, and they are the only terms that survive in the kinetic limit. We define separately the Laplace 
transforms of the irreducible "excitation" diagrams (7Z C ^) and the irreducible "ladder" diagram (£): 



TZT(z) := f dte- tz V A 2W V t (vr 

j~r+ , , r „ 



(5.11) 



M > 2 

7r irr. 



C{z) := f dte~ tz V V t (7r) = / di e-* z V t ({(l, 2)}). (5.12) 

|tt|=1 

Here and in what follows, we omit the specification tt e P under the summation symbol. We observe that, in 
||5.12| |, the only element of V\ is the set containing the single pair (1, 2). The operators TZ^(z) and C(z) have already 
appeared in Theorem l4.4l We will prove Theorem l4.4l in Section l531 First, we establish some helpful estimates 



5.2 Estimates on the Dyson expansion 
5.2.1 A priori estimates 

The following Lemma l5Al is a useful a-priori estimate. Its main assertion, Statement 2), i.e., eq. (|5.14b , gives a bound 
on Vt(7r), the contribution of the irreducible pairing 7r to the dynamics, in terms of the temporal coordinates t. In 
particular, the sum over the other coordinates, x and / is already performed. This is possible because the matrix 
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elements of the free dynamics (e~ ltHs )(0, x) decay exponentially in space, for fixed t; (see Statement 1 of Lemma 
15.11 or eq. (15.171 ). Eq. (15.181 tells us that one can sum over x at the cost of introducing an exponential growth in 
time. This exponential growth in time is also visible in l !5.14b , in the factor e 2 * Ce( - 71 \ However, this exponential 
growth is harmless, because the reservoir correlation functions %/} on the RHS of ( 15.141 are exponentially decaying 
in time, by Assumption [Z2j and the growth constant c e (ji) can be chosen arbitrarily small. In particular, it can be 
chosen smaller than the reservoir decay rate g-R, and this fact will be exploited in Section l5.3.2l . 

Lemma 5.1. Suppose that Assumption \2.1\ holds (with some S e > 0) and define 

Ce(5) := sup fceT d sup| SK |< 5 |3e(fc + k)|, (c e (5) < oo, for < 5 < 5 £ ) 



b d (S) 



-S\x\ 



(b d {6) < oo, for 0<6). 



Then the following statements hold true. 

1) For any k e C d with |Sk| < ji, for some ji < S e , 



|| e i(/» 1 X) e -ite(P) e -i(/»,X)|| < e tc c ( 7 i) ) t > o. 

2) Let 7r e V n , and choose constants < 7 < 71 < 5 e . For any k e C d satisfying \$$k\ £71 — 7, 

M2 7 ) [Hit - 7 - M)] 2 " 2 2n e 2tc ^ 



\\J K vMJ- R \\ < 



/2n-l \ 

/ n d*< n Ws-tr) 

Q=ii<. ..<t 2 n=t \ *= 2 / (r,s)G7T 



We reca/Z f/tat || • || /« ([5.141 refers to the operator norm on 3S{3§2(j%k)). 
Proof. Statement 1) 

Recall that Hg = e(P). By analytic continuation from 3k = to < <5 £ , one has that 

e i( K ,X) e -ite(F) e -i(«,X) _ e -ite(P- K )_ 

Since, for |9f«| < 7x 

|| e - ite (p-«)|| < e t||9fe(p-«)|| < e M-n) j t >0, 
the claim (15.131 is proven. We observe that ([5.131 implies 

\(e- itHs )(x,x')\ <e tc °^e-^ x '- xl , for any < 71 < S e , t > 0, 

and hence 

^ e 7|a: '- a:| |(e- itffs )(a;,a;')| < e tCe(7l) & d ( 7 i - 7), for any < 7 < 71 < <S £ , t > 0. 



(5.13) 



(5.14) 



(5.15) 
(5.16) 

(5.17) 
(5.18) 



Statement 2) 

To estimate the integrand in (15.91 1, we choose < 7' < 71 < 5 e and find that 



E 

y',z' 



>i'{\y'-y\+\ z '- z \) 



^2 CniL x, l_)l X2nMn .. .l X2 . h Ut 2 ^t 1 ^x 1 ,i 1 (y,z;y',z') 



< (sup I X! e 2tCe(7l) (^(7i - 7')) 2 ™, 



(5.19) 



where we can replace "X}/' by 2 2n , the number of terms in the sum. The bound (15.191 1 is obtained by applying 
(15.181 1 2n times. 
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For clarity, we illustrate this with an example: Take n = 4 and (1%, l 2 , 13, 14, 15, Is, h, Is) — {L, R, L, L, R, L, R, R). 
First, we notice that 

\{i Xs j 8 . . .ixa.izUtv-txZxKh) {y, z \y', z ')\ ( 5 - 2 °) 

vanishes unless xx = y and xg — z', and that it is bounded by 

I w(t 3 -ti,x 3 - x x ) x w{t A ~ t a , x 4 - x 3 ) x w(t 6 - U, x 6 - x 4 ) x w(t -t 6 ,y' - x 6 ) 

1 X w(t 2 — 0, X 2 — z) X w(t 5 —t 2 ,X 5 — X 2 ) X w(t 7 — t 5 ,X7 — x 5 ) x w(t s — t 7 , x 8 — x 7 ), 

where w(u,x) := |(e~ itlHs )(0, x)|, t x = 0,t 8 = *• 

We use the decomposition (recall that xx = y and x$ = z) 

W - V\ < \X3 - XI I + |x 4 - X 3 | + |X6 - X4| + |y' - x 6 |, 
|z' - z| < |x 2 - 2; I + |X S - X 2 | + |x 7 - X 5 | + |x 8 - X 7 |, 

and ([5.211 to factorize the sum over y', z' , x on the LHS of l !5.19|l . Those sums can then be carried out with the help 
of j5.18|l , yielding the bound 



/\\8 



M11 - i)) 



exp (c s ( 7l ) [(t a - t 6 ) + (t 8 - ^) + (*4 - *a) + (*s - *i)D 
x exp (c £ ( 7 i) [(t 8 - i 7 ) + (t 7 - t B ) + (t s - t 2 ) + (t 2 - 0)]) 



= (^(7i-7')) 8 e 2tCe(7l) (5.22) 

Note that this bound only depends on \ir\ and t, and not on t, I, or it. Hence it can be applied for all I, which yields 
the factor 2 2 " in (5l9l> . 

For a linear operator VV on l 2 (Z d x Z d ), a straightforward application of the Cauchy-Schwarz inequality yields 

||W|| < b d {25) I sup V |W(2/,z;y',z')|e 5(|y '^ l+|z '" 2|) ■ (5.23) 

Starting from the explicit definition of J K V{(tt) JL re (as in Il3.2b and \5.9\ ) , one uses (|5.23l l and J5.19|l with 7' := 
7 + This yields Statement 2). 

□ 

5.2.2 A combinatorial estimate 

In the next step of our analysis of the Dyson series, we show that one can perform the integration over all pairings 
7r and temporal coordinates t contributing to i5.1H . The following lemma is purely combinatorial, i.e., it only 
employs notions introduced in Definition 15. II 

Lemma 5.2. Consider a positive function h on E + and a pairing it e V. We define 



2n-l 



Xt(7r):= / (n d ^) II h (ts~t r ), withn=\w\ 

0=ti<...<i 2 „=* i=2 (r,s)Eir 



(5.24) 



Then 



XtW < ( £ xt(7r)) x exp ft J + dwh{w)\ , (5.25) 

7rirr. 7rmin. irr. \ K / 

and, if it is the minimally irreducible pairing in V n and z£R, 

f dte~ tz X t{ir) < ( f dwh{w)e~ w A x f f dyf dw h{y + w)e~ w A . (5.26) 
Jm+ \Jr+ J \Jr+ Jr+ J 



22 



Proof. Given tt 6 V n , we can relabel the times ti , . . . , t 2n by setting 

m — t Ti , Vi — t Si for i = 1, . . . , n. (5-27) 
Using our conventions for the labels of the pairs (r^, Si), it follows that 

< Ui < Vi < t, < Ui < Ui + \ <t, = u\,t = max{wi}. (5.28) 

Conversely, a set of n pairs of times i = 1, . . . , n, satisfying (|5.28l l uniquely determines a pairing tt e V n and 

corresponding times = t\ < . . . < t 2n = t. 




t\ t% ta ti ^5 t$ 
U\ u 2 V\ u 3 v 2 v 3 



Figure 5.4: This figure illustrates the change of variables from (n, t), with tt £ V-j, and t± < . . . < ta, to (ui, 1^)^=1,2,3, with Ui < Vi and 

Ui < u i + i. 

Consider an irreducible pairing tt' E Vn' ■ It is easy to see that we can always find a subset j\ , . . . , j n of 
{1, . . . , n'}, for some n < n', such that 

1) the pairs (r^ , s^), i = 1, . . . , n determine a minimally irreducible pairing tt 6 "P„; 

2) these pairs contain the boundary points, i.e. rj 1 = 1 and m&x i {sj i } = 2n'. 

We write tt' — ► 7r whenever 7r and tt' are related in this way; (note, however, that tt is not uniquely determined). It 
follows that 

E xt(^)< E E x*(*0- < 5 - 29 ) 

| tt | — n \tt\ < n 1 7r | — n 

7T irr. 7r min.irr. 7T — > 7r 

For n' > 2, the inequality is strict, since tt is not necessarily uniquely determined by tt', and hence the same irre- 
ducible tt' can appear more than once on the RHS of (|5.29t . 

Using the change of variables 115.271 1, one can convince oneself that, for all m := n' — n > 0, 

„ rn 

E XtW) = Xt{K) j dudv (5-30) 

K'l = o < mi < . . . < u m < t i=1 

where 7r is the minimally irreducible pairing in V n , and where we have abbreviated du := dui . . . du n and du :— 
dvi . . . dv n . The relation l|5.30l l expresses the fact that one can add any set of pairs, corresponding to times u, v 
satisfying the first two conditions of d5.28b , to a minimally irreducible tt, thus obtaining a new irreducible pairing 
(see also Figure 1531 . 

By explicit computation, 

P in „ * „ x rn 

2_\. / dudv_ JJ^ h(vi — Ui) < E] / I / dwh(w) I 



< ui < .. . < u m < t mfc * 0<«i<...<« m <t 

< ni < ti; < t 



< exp ( t / (iw/i(ui) ) , (5.31) 
which proves the bound A5.25I I starting from | |5.29l l and | |5.30t . 
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Ul U 2 Vi V 2 



Ui V\ U 2 V 2 



U\U 2 Vi V 2 



Figure 5.5: Illustration of 15,301 , Three pairings in V5 contributing to the LHS of 15,301 , We have chosen n to be the minimally irreducible 
pairing (1, 3), (2, 5), (4, 6) in V3, as in Figure IB~4l , For each of these 3 pairings in V5, the five pairs («i,«j)i=i,... > 5 contain a subset of three 
pairs identified with n. We have only shown the two other pairs, relabeling them as (ui, 1^)^=1,2- The same strategy is used to prove 15,301 in 
general. 

When we perform the change of variables 1 15.271 1 for a minimally irreducible pairing ir, the variables u, v satisfy 
the constraint ui + i < V, < ui +2 in addition to the constraints < Ui < ui + i < t and < Ui < Vi < t. Let tt be the 
minimally irreducible pairing in V n . Then (ui = is a dummy variable) 

/>oo pvi />00 

dte- tz X t(^)= dvMv! - Ul )e~^ Vl ~ u ^ / du 2 dv 2 . . . 

JO JO Jvi 

CV n — 4 rOO fVn — 2 poo 

du„_ 2 / dv n - 2 ... dii„_! / dv n ^ 1 e^ z(v "- 1 ^ v "- 2 ' ) h(v n ^ 1 - ii„_i) 

'fn-5 Jv n -3 J V n -2 



Vn-1 



du n dv n e- z ^- v ^h(v n -u n ). (5.32) 



Performing the change of variables Wi = Vi — Vi-i and y,; = — Ui (for i > 1) and extending the range of 
integration of yi to R, the above expression factorizes and one obtains the bound (|5.26t . 

□ 

5.3 Proof of Theorem IOI 

In this section, we prove Theorem l4.4l Statement 2) is proven separately for C(z) and lZ c ^(z) in Sections 15.3.11 and 
15.3.21 respectively. Statement 3) is proven in Section 15 .3 . 1 1 and Statement 1) in Section l5.3.3l 

It is mainly in Section 15.3.21 that we use the preparatory work summarized in Lemma 15.11 and Lemma 15.21 in 
order to obtain a bound on TZ c ^{z). 

5.3.1 Properties of C{z) 

We compute C(z) starting from 1)5. 9b and 1 15.121) : 



C(z)= f dte~ tz V V l v ,JA t li, x 

JWL+ , ,,,-rr „1 



$(t) I = L 

(5.33) 



xez d i,i<£{L,R} { ip{-t) l = R 

To display the result, we introduce the functions ip+ , ip- as 

tp+(z)= [ dt^{t)e itz , i>-(z)= I dt^{t)e itz , z e C, (5.34) 
Jr+ 



with ip as defined in Section l2.3.21 (we recall that tp(u) decays exponentially). Since ip{—u) = ip(u) (as follows from 
42.181 1), one has that 

tl>+(z) = i(>-(z), ^(«)=V+(«)+^-(*). with|3z|< 5R . (5.35) 
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Using d533l . we calculate C{z)S, for S e L 2 (T d x T d ), 



0C(*)S)(* + §,*-§) 



= dk' (i, + [e{k - |) - + |) + iz] + V- + f ) - <tf - §) - N) + | , fc' - f ) ( 5 .36) 
^ dfc' (V+[e(fc - |) - e(fc' + |) + iz)] + ^_[ £ (fc + |) - e(fc' - |) - iz]) S{k + |, fc - |). 



The claim about £(z) in Statement 2) of Theorem 14.41 follows by noticing that the above expression can be 
analytically continued in z and p. This follows from the analyticity of e ( Assumption ^ . 1 b and ip+ , ip- (consequences 
of Assumption l2.2b . 

To prove Statement 3), we first check that (£(0))o = M° by setting p = and z = in (|5.36l l, and using Il5.35b . It 
remains to verify that 

X 2 (M K - M°) = iA 2 (k, Ve) = (ad(iff s )) A 2 K + 0((\ 2 k) 2 ) (5.37) 

as operators on L 2 (T d ), where (k, Ve) is the multiplication operator given by the function k i— > (k, Ve)(fc). Equation 
l!5.37b follows by writing explicitly 

((ad(iffs))p«)(A)=i(e(A+|)-e(*-|))e(A), 6eL 2 (T d ,dk), (5.38) 
expanding in powers of p and putting p = \ 2 n. 

5.3.2 Properties of TZ c x x (z) 

Choose positive constants 71 > 7 > 0, as in Lemma [5. 1[ and define the quantity Xt{^) as in Lemma [5. 2[ with h 
given by 

h(t) :=2 2 6 d (7 1 -7-|3 K |) 2 A 2 |^(t)|. (5.39) 
It follows from Statement 2) of Lemma IBTTI and equations i5.9) , l !5.11|l that 

\\J K 1lf(z) J- K \\ < 6 d (2 7 ) / die**^-*** ( V xtW). (5.40) 

7r+ V W>i / 



Ki > 

7r irr . 



and hence, using Lemma [5721 that 

||J K ^r(^)J_ K || < 6 d (2 7 ) / dte-C**-* V XtW, with a:=2c £ ( 7l )+/ dw/iM 



7r min. irr. 



/R+ 

< b d (2j)( f dw h(w)e~ w ^ z ~ a A F ( [ dy f dw h(y + w)e- w ^- a ) 
\Jm+ J \Jr+ Jr+ 

where F{x) := j^— , provided that |sc| < 1. To prove the first inequality above, we use (|5.40b and J5.25|l , and, for the 
second inequality, we use ( 15.26b and sum the geometric series. 

Statement 2) of Theorem 14.41 now follows by fixing the constants and using the exponential decay of ip. For 
example, choose 71 , 7 such that 

2c e (7i) < \qr, 1 ■= \li (5.41) 

and 8' 2 small enough such that for | A | < 8' 2 

[ dwh(w)<^-g K , [ dy [ dw h(y + w)e- w ^i gR ^ < 1. (5.42) 
Jr+ 4 J R + J R + 

Then J4.47|l is satisfied with 6[ := 771, g' := j9r an d 8' 2 as determined above. 
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5.3.3 Proof of equation J4.44b in Statement 1) of TheoremlPl 

To simplify the following calculations, we abbreviate 

TZT(z) := TZT(z) + \ 2 C{z), K s (z) := (z - ad(Lffs))- 1 . (5.43) 

By the self-adjointness of ad(-ffg)/ ° ne nas that ||7£s(z)|| < |3fJz| _1 . We choose A and z such that Sftz > and 
\\TZY(z)TZ s (z)\\ < ISRzr 1 !!^"^)!! < 1. Then 

-tz syX 



K x (z) := / dte~ tz Z, 



= n s (z)(KY(z)n s (z)y 



nGZ+ 



= n s (z)(i^'R i r(z)n s (z)y 1 

= (z-ad(ii7 s )-^rW) _1 

= (z-ad(ii7 s )-A 2 £(z)-^rW) _1 , (5.44) 

where the second equality follows by Laplace transforming (|5.101 , and the third equality represents the sum of a 
geometric series. Hence, Statement 1) of Theorem l4.4l is proven. 



6 Proof of Theorem 113 



In this Section we prove Theorem 14.51 Our reasoning is based on a standard application of analytic perturbation 
theory and the inverse Laplace transform. 
We abbreviate 

A(z, A, k) := (ad(Lffs) + A 2 £(z) + K cx (z)) a2r - \ 2 M K (6.1) 

and we define 

G := { (z, A, k) e C x R x C d Kz > -g', \k\ < <5 ki „, |A| < win(S' 2 , J A-) \ , (6.2) 
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kin 



where g' , <5 l7 are as described in Theorem 14.41 and S^ n is as described in Theorem l4.3l Theorem 14 .41 implies that, 
on the domain G, the function \ 2 M K + A(z, A, k) is analytic in the variables (z, n) and, for 5Rz large enough, 

{Kx(z)) X 2 K = (z- \ 2 M K - A(z, A, k))- 1 . (6.3) 

We may extend the (operator-valued) function z i— » (R-\(z))\2 K into the region 5Rz > —g'. This will be useful, 
because, at the end of this section, we calculate the reduced evolution (Z ( A ) A 2 K from the inverse Laplace transform 
of (TZ\(z)) X 2 K . From J6.3I I we see that any singular point of the function z i-» (lZ\(z)) X 2 K must satisfy 

z € sp(X 2 M K + A{z,X,k)). (6.4) 

Recall that by Theorem 14. 31 M K has a simple isolated eigenvalue /ki n (t), and let f2 c C be as defined in (|4.30t , i.e., 

Q:= U (s P Af K \ {/ki„(K)}) ■ (6.5) 

N<<5 ki „ 



The following two lemmas describe the singularities of (Z t A 



Lemma 6.1. There is a constant c\ and a function c(A) with c(A) \ 0, as A \ 0, smc/? that, for any z satisfying J6.4| >, one of 
the following two statements holds 

dist(z,A 2 0) < A 2 c(A), or dist(z, A 2 / kin (K)) < ciA 4 . (6.6) 
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Proof. From Theorem 14. 41 we infer that 

\\A(z,X,k)\\ = A 2 ||(£(z)-£(0)) A 2 K || +0(A 4 ) + 0((A 2 k) 2 ) asA\0,A 2 K \0, (6.7) 

with (C(z) — C(0))\2 K bounded and analytic in (z, n) on G. Since M K is bounded, there is a constant r(m) > 0, for 
all to > 0, such that 

sup -M*)- 1 1| < to. (6.8) 

z£C, dist(z,spM K )>r(m) 

Choose to -1 := sup( z A ^ eG A~ 2 ||v4(z, A, k)|| (by l|6.7|l , mT 1 = 0(A°)). Using the Neumann series for (z — X 2 M K — 
A(z, A, k))~ x , it follows that, if dist(z, A 2 spM K ) > A 2 r(m), then z cannot satisfy 116.411 . 

If, however, dist(z, A 2 spM rt ) < \ 2 r(m), then \\A(z, A, k)\\ = 0(A 4 ), as A \ 0; (this follows from $67$ and the 
analyticity of C(z)). The claim now follows from analytic perturbation theory, using that A 2 /ki n («) is an isolated 
simple eigenvalue. □ 

Lemma 6.2. For sufficiently small | A|, there is a unique z —: z at a distance O(\ 4 )from A 2 /ki n («) satisfying j6Aj . Let P X > K 
be the residue of [z — A 2 M K — A(z, A, at z — z. It follows that P A,ft is a rank one-operator and 

\\P x > K -P£J = 0(\ 2 ) (6.9) 

with Pg in the one-dimensional spectral projection ofM K corresponding to the isolated simple eigenvalue /kin («), as in Theorem 
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Proof. By analytic perturbation theory, the operator \ 2 M K + A(z, A, k) has at most one eigenvalue at a distance 
0(A 4 ) of /kin(^)- This means that 116.41 has at most one solution at a distance 0(A 4 ) of /kin(«0- We now prove that 
there is at least one solution. Indeed, if no such solution existed, we could choose a contour 

C K , a = {zeC\\z- f kin ( K )\=a}, a>0, (6.10) 

with a small enough such that C K . a stays away from Q. We then calculate 

27ri(P k K in - 0) = / dz(z — A 2 Af K )~ 4 — f dz(z — \ 2 M K — A(z, A, k)) _1 

J\ 2 C K . a J\ 2 C K . a 

dz(z - X 2 M K )~ 1 (1 - (1 - A(z, A, k){z - A 2 M K )- 1 )- 1 ) 

A 2 C K , a 

where 

b(a,K) := sup \\(z — M re ) _1 ||, 

26C«, B 

and, here and in what follows, the contour integrals are meant to be oriented clockwise. Since the last line of 1)6.111 
is of order A 2 , we arrive at a contradiction to the fact that P£ in ^ 0. 

The claim about the residue is most easily seen in an abstract setting: Let F(z) be a Banach-space valued analytic 
function in some open domain containing 0, and such that € spF(0) is an isolated eigenvalue. We have hence the 
Taylor expansion 

F(z) = y2-rF n , F n :=F<™)(0), e s P F (6.12) 
* — ' n! 

If || F\ — 1 1| is small enough, then also Ff Fo has as an isolated eigenvalue. We denote the corresponding spectral 
projection by lo(Ff" Fq) and we calculate 

Res(F^)- 1 ) = Res(F + zFx)- 1 = (Res(Ff X F + z)- 1 ) Ff 1 = l (Ff 1 F )Ff 1 . (6.13) 

The last expression is clearly a rank-one operator. In the case at hand, Ff 1 = 1 + 0(A 2 ), as A \ 0, which yields 
(62). □ 
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Figure 6.1: The (rotated) complex plane. The black dots and thick black line indicate the spectrum of A 2 M°: The upper dot is the eigenvalue 
and the thick vertical line is the continuous spectrum. In the picture, we have drawn only one other eigenvalue, but, in general, there can be 
more than one (or none) further eigenvalues. The function X 2 M K + A(z, A, k) is analytic above the lowest gray (rectangular) region. The other 
gray regions contain the singularities of the function (7l\(z)) X 2 K for (z, X, re) £ G. The integration contours T, T' and A 2 C K , a , A 2 C are drawn 
in dashed lines. In this picture, the contour A 2 C K . a encircles A 2 /(A, re), for all (A, re), (i.e., such that (z, A, re) £ G), which can be achieved by 
choosing a large enough. 



We set /(A, k) := z and we define P X,K as the residue of (z — \ 2 M K — A(z, A, at z = z. It is clear that /(A, k) 

and P x - K enjoy the analyticity properties claimed in Theorem l4.5l 
We define the horizontal contours 

r := {z e C\z = 1 + iR}, V :={ze C\z = ~(g' - e)+iM}, (6.14) 

with I large enough such that all singular points of z i— » (lZ\(z))\2 K lie below V, and e > small enough such that 
all singular points with 5Rz > — g' lie above T' (the notions 'below' and 'above' are meant as in Figure fSTTt . These 
contours are oriented from left to right. By Theorem l4.31 we can construct a contour C which encircles 51 and such 
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that /kin(«0 is separated by a gap g from this contour: 

g := inf 5R/ kin ( K ) - sup^C > 0. (6.15) 

|«|<<5ktn 

By performing an inverse Laplace transform we find that 

(Z ( a ) a2k = ^ J dze t *(z-X 2 M K -A(z,X, K ))- 1 . (6.16) 
r 

For A small enough, Lemma [6.11 ensures that one can deform contours and obtain 

1= J + J + J- <617) 

r A 2 c a , K a 2 c r' 

The first term on the RHS of ([6.17l l equals e tA2 /( K A )p A «; this follows from Lemma [6721 The second term is domi- 
nated by 

e A 2 t(s upGRc')) J ^\\^-x 2 M k )- 1 \\\\{1-{1-A{z,\k){z-X 2 M k )- 1 Y X \V (6-18) 

A 2 C 

By the choice of C' A and the bound i6.7h the integral on the RHS is bounded by a constant, for A small enough. 
The third term of the RHS of ([67171 is split as 



dz e (z — \ M K — A(z, A, k)) = / dze (z — X M K )~ (6.19) 

r' 

+ / dzc tz (z- X 2 M k )~ 1 A(z,X,k)(z- X 2 M k - A(z,X,k))- x . 



The first integral can be closed in the lower half -plane and equals 0, the second integral has an integrand of order 
z~ 2 for large z, and hence its contribution is bounded by a constant times er t{ - 9 ~ e \ 
It follows that the crucial estimate 04.491 1 holds with 8\ := <5km an d g as in ([6.151 1. 



APPENDIX A 

Here we consider the effective structure factor, which, in Section [727731 has been introduced as the Fourier transform 
of the reservoir correlation function. 

We use the spectral theorem to represent the positive operator uj as multiplication by £ G M + . There are Hilbert 
spaces f)^ for £ e M + such that f) = f mu+ d£f)£, and for all ip 6 F), there are <p£ e f)^ such that 



¥> = / u)ip= I d^ip^. (A-l) 

The structure factor G f) has been introduced in Section [727751 We construct an effective form factor i/)' 3 as an 
element of t) © (). We choose F)_£ to be isomorphic to (){, and we define = f ffiR 0^ as an element of f) © f) ~ / R f)£ 
by setting 



A /3 ._ J V 7 ^ 5 -! 



e>o, 



(A-2) 

The function 0^ plays the role of the form factor if one constructs the positive-temperature dynamical system. We 
just note that 

m = \\4\\li- ( A - 3 ) 
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Assume that the on-site one-particle space is given by f) = L 2 (M. d ), and the one-particle Hamiltonian acts by 
multiplication with a function £(r), where r := \q\, for q E R d . We also assume that r i— ► £(r) is differentiable and 
monotonically increasing. Hence we can define the inverse function £ i— > r(£). The form factor € L 2 (K d ) is taken 
to be spherically symmetric, </>(g) = 0(r). Then the Hilbert spaces f)£ are naturally identified with i 2 (S d_1 ), and 

5 / \ (l-e«)-'/^( r (- f )), £<0, 

where lgd-i € L 2 (S d_1 ) is the constant function on § d_1 with ||lgd-i || = 1. 

Next, we return to Assumption ^. 21 By properties of the Fourier transform, e.g. Th. IX.14 of [19], this assumption 
is equivalent to the assumption that ip extends to an analytic function in the strip |3ff | < <?r, and 



sup / dx \tf)(x + iy)\ < oo. (A-5) 

-gn<V<9n Js. 

Starting from expression dA-4b , one can check condition jA-5\ in concrete examples. E.g., for a relativistic dispersion 
law, £(r) = r, \A-5\ is satisfied whenever 

sup / dx \x + iy\ d ~ 2 \(/>(x + iy)\ 2 < oo. (A-6) 

-9r<w<sr Jr 
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